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We introduce a model of a dual-core optical waveguide with opposite signs of the group-velocity- 
dispersion (GVD) in the two cores, and a phase-velocity mismatch between them. The coupler is 
embedded into an active host medium, which provides for the linear coupling of a gain-loss type 
between the two cores. The same system can be derived, without phenomenological assumptions, 
by considering the three-wave propagation in a medium with the quadratic nonlinearity, provided 
that the depletion of the second-harmonic pump is negligible. This linear system offers an optical 
realization of the charge-parity (CP) symmetry, while the addition of the intra-core cubic nonlin¬ 
earity breaks the symmetry. By means of direct simulations and analytical approximations, it is 
demonstrated that the linear system generates expanding Gaussian states, while the nonlinear one 
gives rise to broad oscillating solitons, as well as a general family of stable stationary gap solitons. 
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I. INTRODUCTION 

Charge-parity-time (CVT) symmetry is the most fundamental type of symmetry in quantum field theory [2, '5], 
where it holds for all relativistically invariant systems obeying the causality principle. Its reduced form, viz., the CV 
symmetry, is almost exact too, save the small violation by weak nuclear forces [4j. The CVT operator is composed 
of three factors : parity transformation, V, which reverses the coordinate axes; charge conjugation, C, which swaps 
particles and antiparticles; and time reversal, T. 

The proof of the CVT and CV symmetries (when the latter is relevant) applies to Hermitian Hamiltonians ( H ), 
subject to the condition H = H', which guarantees that the spectrum of the Hamiltonian is real. However, one cannot 
deduce from the CVT or CV symmetry that the respective Hamiltonian is necessarily Hermitian [5]. Indeed, the 
consideration of Hamiltonians which commute with a reduced symmetry operator, VT , demonstrates that they may 
contain an anti-Hermitian (dissipative) part, provided that it is spatially antisymmetric (odd), while the Hermitian 
one is even [6]. The spectrum of such a Hamiltonian remains purely real up to a critical value of the strength of the 
anti-Hermitian part, at which the VT symmetry is broken, making the system an essentially dissipative one (recently, 
a model with unbreakable VT symmetry was found; it includes defocusing cubic nonlinearity with the local strength 
growing from the center to periphery [56]). 

While in the quantum theory the possibility of the existence of non-Hermitian VT- symmetric Hamiltonians is 
a purely theoretical one, such systems have been realized, theoretically [M3 and experimentally [Hl - BD] . in optics, 
making use of the fact that the wave-propagation equation, derived in the standard paraxial approximation, is identical 
to the Schrodinger equation in nonrelativistic quantum mechanics. In this context, the spatially symmetric and 
antisymmetric Hermitian and anti-Hermitian terms of the Hamiltonian are represented, respectively, by even and odd 
distributions of the refractive index, and of the local gain-loss coefficient in the photonic medium. A VT- symmetric 
electronic circuit was built too, following similar principles }21j . 

The essential role played by the Kerr nonlinearity in optics has suggested the development of models in which the 
Hamiltonian includes a quartic Hermitian part too. The nonlinearity gives rise to families of VT -symmetric solitons, 
that were investigated in detail in continuous and discrete systems including VT -symmetric dual-core 

couplers [281130] . Models combining the VT symmetry with quadratic nonlinearity in the dynamical equations (i.e., 
cubic terms in the respective Hamiltonians) were also elaborated [SIMM] . 

As a subject of quantum field theory, the CVT and CV symmetries mainly relate to elementary particles [MUST] . 
On the other hand, the above-mentioned works on the implementation of non-Hermitian 'P'T-symmetric Hamiltonians 
in photonics suggest looking for a possibility to design optical settings that would realize non-Hermitian Hamiltonians 
featuring the full CVT symmetry, as well as its CV reduction. A possibility to implement the former symmetry was 
recently explored in Ref. [58]. which addressed not optics, but rather a two-component Bose-Einstein condensate 
with the spin-orbit coupling between the components, one of which is subject to the action of loss, and the other 
one is supported by gain. In terms of optics systems, the symmetry of that models is similar to the VT symmetry 
in a dual-core waveguide, with a combination of continuous V transformation acting in the longitudinal direction, 
and another V transformation which swaps the two cores. A similar symmetry was proposed in Ref. [30] , which put 
forward a "PT-symmetric coupler subject to the action of “management”, in the form of periodic simultaneous switch 
of the signs of the coupling and gain-loss coefficients. 
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The present work aims to offer emulation of the CV symmetry in a two-component optical system, which, at 
the phenomenological level, may be considered as a dual-core waveguide with opposite signs of the group-velocity 
dispersion (GVD) in the cores and a phase-velocity mismatch between them, embedded into an active medium. We 
demonstrate that the system can be derived, without phenomenological assumptions, as a model of the spatial-domain 
propagation for two fundamental-frequency (FF) components with orthogonal polarizations of light, pumped by an 
undepleted second-harmonic (SH) wave in a birefringent medium with the \^ nonlinearity. We further investigate 
conditions for persistence and breaking of the CV symmetry, both analytically and numerically. In particular, the 
addition of cubic (Kerr, alias x®) terms to the system with the active coupling breaks the symmetry at the nonlinear 
level, but helps to stabilize confined breather states (oscillatory solitons), and gives rise to a family of stable stationary 
gap solitons. 

The paper is organized as follows: Section 2 introduces the model, in its both forms (phenomenological and the one 
based on the x^ interaction) and reports analytical results. Numerical findings for the linear and nonlinear systems 
are presented in Section 3, and Section 4 concludes the paper. 


II. THE MODEL AND ANALYTICAL RESULTS 


A. The system: phenomenological formulation 


At the phenomenological level, we consider the copropagation of optical modes u and v in a dual-core coupler with 
opposite GVD signs in the cores (cf. Refs. UliTl - fTTI . where a similar feature was introduced in different contexts, and 
also Refs. pT?lT]~T| . where systems with opposite signs of group velocities were considered in the contexts of coupled 
right- and left-handed waveguides), and a phase-velocity mismatch 2 q between them. The linear coupling between 
the cores is provided by cross-gain terms, with strength 7 , which is possible when the coupler is embedded into an 
active medium, as recently proposed in Ref. [35] . The model is represented by the following system of propagation 
equations, which include the Kerr nonlinearity too, with respective coefficient a (all the quantities are dimensionless): 

iu z + (l/2)u t t - qu + cr\u\ 2 u = i'yv, ( 1 ) 

iv z - (l/2)v tt + qv + a\v\ 2 v = iju. ( 2 ) 


Here 2 is the propagation distance, t is the reduced time [37] > the GVD coefficients are scaled to be ±1, and 7 < 0 may 
be transformed to 7 > 0 by changing v —> —v. Positive and negative values of a can also be transformed into each 
other by substitution (u*,v*) = ( v , u), therefore, in what follows below we consider only a > 0. Then, rescaling allows 
one to fix a = 1 , but we prefer to keep it as a free parameter, the variation of which helps to monitor a transition 
from the weakly nonlinear system to a strongly nonlinear one. 

It is relevant to mention that a dissipative discrete system with opposite signs of the discrete dispersion and a 
wavenumber mismatch between the components was introduced in Ref. |46| . However, that model included the 
dissipative coefficient in a single equation, therefore it did not realize the symmetry considered here. 

Equations |l]) and © can be derived from the non-Hermitian (complex) Lagrangian, which is usual for VT- 
symmetric systems: 


L = 



i (■ u*u z + v*v z ) + ^ (ju t | - |u t | 2 ) + q (H 2 

/ +00 

(u*v + uv *) dt , 

-OO 



(3) 


which generates the respective non-Hermitian Hamiltonian in an obvious way. The total energy, 

/ +00 

[| u(z,t)) I 2 + I v(z,t) I 2 ] dt = E u (z) + E v (z), 

-OO 

is not conserved by Eqs. Q and ([ 2 ]). Instead, the system gives rise to the following energy-balance equations: 

/ -(-00 

(■ uv * + u* v) dt. 

-OO 


d.E u dE v _ 1 dE 
dz dz 2 dz 


(4) 


(5) 


The equality of dE u /dz and dE v /dz, i.e., the conservation of E u — E v . means that the linear coupling of the present 
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type causes mutual amplification or attenuation of both components. 

The linear version of Eqs. 0 and ([2]), with a = 0, are invariant with respect to the CV transformation, defined as 

( u , v) (u = v* ,v = u *), ( 6 ) 

where the swap of u and v stands for V, and the complex conjugation - for C (conserved E u — E v may be considered 
as the respective charge). It is relevant to compare the system of Eqs. 0>0 and their invariance transformation 0 
with the previously studied model of the "PT-symmetric coupler, which was based on the following equations (28U30[ : 

iu z + (l/2)u tt + cr\u\ 2 u = i'jv, (7) 

iv z + ( 1/2 )v t t + cr\v\ 2 v = —vyv. ( 8 ) 

Obviously, Eqs. ([7]) and 0 are invariant with respect to transformation ( u,v,z ) —> (u = v*,v = u*,z = —z), which, 
in the present context, may be considered as corresponding to the CVT symmetry, the reversal of z playing the role 
of additional T. 

It is relevant too to compare the present model to the system of equations with opposite GVD terms, coupled by 
the usual conservative terms, rather than by those representing the gain and loss m- 

iu z + (\/2)utt + a\u\ 2 u + Kv = 0, 

iv z - {l/2)v tt + a\v\ 2 v + Ku = 0, ' ' 

where K is a real coupling constant. The linear version of this system is invariant with respect to the anti-CV 
transformation: (u,v) —► (u = v* ,v = — u*) , “anti” corresponding to the relative sign flip, cf. Eq. 0. 

The nonlinearity breaks the symmetry of system 0 , §. as the opposite relative signs of the GVD and cubic 
terms in the two equations make it impossible to swap u and v, which represents the V transformation in Eq. §. 
Nevertheless, nonlinear effects are obviously interesting too. It is demonstrated below that the nonlinearity creates 
solitons in the present system. In this connection, it is relevant to mention recently introduced nonlinear models 
with alternating gain and loss, which do not obey the condition of the VT symmetry, but nevertheless support stable 
solitons HEHSill- 

A solution to the linear version of Eqs. 00 in the form of plane waves, {u, u} = {uo, r'o} exp (ikz — icot ), produces 
a dispersion relation for the wavenumber and frequency: 




( 10 ) 


Obviously, in the case of q > 0 the spectrum given by Eq. (10) is pure real, provided that 


< 7thr = q, 


( 11 ) 


while in the case of q < 0 the spectrum always includes an imaginary component. The change of the spectrum from 
real to a partly imaginary one, with the increase of the gain-loss coefficient, at | < y| = q (provided that q > 0 ) implies 
the breakup of the CV symmetry, similar to the phase transition which is the generic feature of Vl ~-symmetric systems 
m- If condition (|II| holds, the spectrum given by Eq. (|10| features a bandgap , 


7 2 - 2 2 

k < q - 7 . 


( 12 ) 


B. The linear model: physical derivation 

While the system of Eqs. 0, § was introduced above phenomenologically, its linear version can be derived, in 
the spatial domain (rather than in the temporal one), starting from the fundamental propagation model for two FF 
and one SH components of light waves, u, v and w, respectively, in the dissipation-free medium with the Type-II x^ 
interaction [5^1 - 155j : 


iu z + ( 1 / 2 )u xx - qu = -v*w , 
iv z + (l/2)v xx - qv = -u*w, 
2 iw z + (l/2)w xx = - ( 1 / 2 ) uv, 


(13) 

(14) 

(15) 
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where x is the transverse coordinate, and q is an FF-SH wavenumber mismatch. Then, adopting the usual approx¬ 
imation for parametric down-conversion, of undepleted SH pump, we replace it by a constant, w = —i'y, neglecting 
Eq. (15), denote v* = v. and apply the complex conjugation to Eq. (14): 


iu, + ( 1 / 2 )u xx - qu = i'yv, 
iv z - (l/2)v xx + qv = iju. 


(16) 

(17) 


These equations differ from the linear version of Eqs. ( IT]). 

As for cubic terms, they can be added to Eqs. (13] 15) as ones accounting for the Kerr nonlinearity in the \^ 
waveguide. However, in terms of Eqs. and © the resulting cubic terms will be different from those adopted in 
Eqs. 0 and 0 : as the above-mentioned complex conjugation of Eq. ( |14| ) will produce the cubic term in Eq. © 


([ 2 ]) only by the replacement of t by x. 


with the sign opposite to that in Eq. ([ 2 ]) [incidentally, in this case the cubic terms do not break the CV invariance of 
Eqs. ([T]) and §]• Furthermore, because the original components, u and v, correspond to two orthogonal polarizations 
of the FF wave, the nonlinear extension of Eqs. ( p~ 6 ] ) and © should also include the respective XPM (cross-phase- 
modulation) terms, viz., (2/3)ajc| 2 u and —(2/3)er| u\ 2 v, respectively, assuming that the four-wave mixing terms may 
be neglected, as usual, due to sufficiently strong birefringence m ■ In the present work, we focus on the nonlinear 
terms adopted in Eqs. 0, 0, while those corresponding to the derivation for the x ^ 2 system will be considered 
elsewhere. 


C. The analytical approximation for broad pulses 


The system based on Eqs. 0,0 can be investigated in an analytical form for broad small-amplitude pulses, with 
widths (t) and amplitudes satisfying conditions 


r 2 » 1/q ; Ug, Vq < 1/ (erg). 


(18) 


In this case, the linearized version of the system yields, in the lowest approximation, two different relations between 
the field components: one solution has 


and another one features 


v(z,t) = 


u(z,t) = 


27 


' + vV - 7 2 


u(z,t ), u(z,t) = e 


= p-V? 2 -') 


—Z7 

q + y/q 2 - 7 2 


v{z,t), v(z,t) = e 


= 




u{z,t), 


v{z,t), 


(19) 


( 20 ) 


where u(z,t) and v(z,t ) are slowly varying amplitudes, in comparison with exp (^±i\/q 2 — 7 2 z'j. The substitution 

of these expressions into Lagrangian (|3]) leads to a real effective Lagrangians for the slowly varying functions (its 
imaginary part cancels out in the present approximation), which, in turn, give rise to either one of the two following 
nonlinear Schrodinger (NLS) equations for the slow evolution: 


d 


± l n — 

± 2 ° eS dt 2 \ v 


^eff 


|u | 2 u 

\v \ 2 V 


= 0 , 


dz \ v 

where + and — pertain to u and v , respectively, while the effective GVD and nonlinearity coefficients are 


D e fi = 


Vq 2 - 


T 


q 


^eff — & 


2q 2 - q 2 


(9 + \7 


2 _ ^2 


( 21 ) 


( 22 ) 


Note that Eq. (21) implies that the system conserves the total energy in the present approximation, which complies 
with the fact that the substitution of relations (19) and (20) into energy-balance equations ([5]) yields dE u ^ v /dz = 0. 
Fundamental solutions to the linear version of Eq. (21) are well known in the form of expanding Gaussians (coherent 
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states, in terms of quantum mechanics) |47| : 




( 23 ) 


where to is an initial width, and ( Uq,Vq ) are arbitrary amplitudes. This result is drastically different from that 
obtained for broad pulses in the linear version of Eq. § with the conservative coupling HU- The latter system gives 
rise to an effective equation for a slowly varying function in the form of a single linear Schrodinger equation with 
the GVD term periodically (in z) changing its sign, thus generating robust oscillating Gaussian pulses, rather than 
expanding ones (23). 


Further, the full nonlinear equations (21) for u and v give rise to commonly known solutions for bright and dark 


solitons, respectively. In particular, the bright NLS solitons with an arbitrary (small) propagation constant, 0<K<g, 


are 


Usd = ■\/2«/c r effe t ' tz sech (VWArfU 


(24) 


It is also worthy to note that the equations for u and v are, severally, Galilean invariant ones, i.e., the linear and 
nonlinear versions of the equation for u give rise, severally, to moving Gaussians and bright solitons, while the full 
underlying system 0 >® does not feature the Galilean invariance. 


D. Gap solitons 


Inside of bandgap (12), it is natural to look for stationary gap-soliton solutions [3?| of Eqs. ([I]), <[2j) in the form of 

{u,v}=e lkz {U(t),V(t)} 1 (25) 


for which Eqs. 0 and 0 reduce to a system of stationary equations, 

~kU+ ~qU+ \U\ 2 U = ijV, 

~ kV ~l d W +qV+ ^ V = ilU 


(26) 

(27) 


(recall a = +1 is fixed). Unlike the usual coupled-mode system for Bragg gratings [57], there is no substitution which 
could reduce Eqs. (26) and (27) to a single equation, therefore the system should be solved numerically, in the general 


case. The stability of the gap solitons should then be tested numerically too. As concerns broad solitons (24), they 
actually correspond to the gap solitons at values of k close to the bottom of bandgap ( 12 ), k = — sjq 1 — 7 2 + k, cf. a 
similar relation between the general gap solitons and broad ones in the standard model of nonlinear Bragg gratings 


A gap-soliton solution to Eqs. (26) and \27\ with strong asymmetry between the two components can be found in 
an approximate form for the limit case of weak gain-loss coupling, y 2 <C q 2 , at k = 0 , i.e., exactly at the central point 


of bandgap (12). The zero-order approximation (for 7 = 0) is 


U = (q + k )sech ^\/2 (q + k)t^j , V = 0, 


(28) 


for q + k > 0 (in this approximation, k = 0 is not required). Then, the first-order correction is determined by the 
linearized equation for V : 


(q -k)V- \ = ilV 2 {q + k )sech (\/2 (q + k)tj . 


(29) 


An exact closed-form solution to Eq. (29) can be found, by means of the Fourier transform, at the center of the 
bandgap, i.e., at k = 0 (a similar solution was reported, in another context, in Ref. my- 

V(t) = i'Yy/^q |y/2qt exp \/2qtj + cosh (^^/2qtj In 1 + exp^—2y/2 qtj j. 


(30) 
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An exact solution to Eq. (29) is available too in the special case of k = —(3/5 )q < 0: 


v{t) = j exp +exp 


— exp ( 8 \/ —t ) arctan 


exp ( -2 A /-t 


— arctan 


exp ( 2 A /-t 


(31) 


These solutions exist under exactly the same condition, q > 0, which was adopted above. Although it may not be 
immediately obvious, both solutions (30) and (31) are even functions of f, exponentially decaying at |<| —>• oo. These 
approximate analytical solutions are compared with their numerically found counterparts below, see Figs. E3 and [16] 


III. NUMERICAL RESULTS 


A. The linear system 


Equations (jlj|2j) were solved numerically by means of the split-step Fourier-transform method [47 
sets of initial conditions. Two of them were taken with a Gaussian pulse in either component: 

u(z = 0, t) = exp(—0.05f 2 /2), v(z = 0,t)=0, 
u{z = 0, t) = 0, v(z = 0, t) = exp(—0.05f 2 /2), 


for four different 

(32) 

(33) 


and two other initial sets are given by Eqs. ( 19p0 ) and (23) with z = 0. The interval for the temporal variable 
was fixed as —800 < t < 800, to ensure that reflections from its boundaries did not affect the results, moreover all 
quantities plotted are dimensionless. 

We start the analysis with the case of weak coupling, 0 < 7 < g e 1. A typical numerical solution for this case, 
presented in Fig. |TJ shows the expansion of the Gaussian launched in the form of initial conditions (32). The v 
component remains weak as the coupling constant is small and, accordingly, the total energy remains very close to 
the initial value. A detailed comparison with Eqs. (19][2T) demonstrates that the asymptotic stage of the evolution, 
for broad pulses, is accurately predicted by the analytical approximation. 

Next, we consider the situation close to the CPT-symmetry-breaking threshold 0 , namely, with 7 = 0.9 for q = 1. 
The respective numerical solution, generated by initial conditions (32), is displayed in Fig. [2j which, naturally, demon¬ 
strates strong coupling between the two components and more dramatic evolution. In this case too, the asymptotic 
stage of the evolution for broad pulses is correctly predicted by the above-mentioned analytical approximation. 

To test the symmetry of the system, we have also performed simulations of the evolution starting from initial 
conditions (33), with swapped components u and v. Comparison of the respective results, displayed in Figs. [3](a) and 
§b), with their counterparts shown above in Fig. [2] confirms the symmetry. Furthermore, the detailed comparison 
of the real and imaginary parts of the two components in both cases (not shown here in detail) exactly corroborates 
the full CV symmetry implied by definition ([6]). 

For 7 > q, when the the CV symmetry of the system is broken, according to Eq. 0, direct simulations (not shown 
here) demonstrate blowup of solutions, as should be expected above the symmetry-breaking point pjlISTj. 

The analytical approximation for broad pulses, based on Eqs. @ and ( [23| , was directly tested by comparing its 
predictions with the numerically simulated evolution commencing from the initial conditions produced by Eqs. (19) 
and (23) with z = 0 and the upper sign in the latter equation. Figure [4] shows that the respective analytical and the 
numerical results are almost identical. The comparison produces equally good results (not shown here in detail) if 
the initial conditions are taken, instead, as per Eq. (20) and Eq. ([23]) with the lower sign, at z = 0. 


On the other hand, for strong coupling, e.g., for 7 = 0.9, when exp ^ ±i\/q 2 — 7 2 z^j is no longer a rapidly oscillating 
carrier in comparison with slowly varying u and v 1 see Eqs. (19) and (20), the analytical approximation is no longer 


relevant. The comparison with the numerical results corroborates this expectation (not shown here in detail either). 


B. The nonlinear system 


Simulations of the nonlinear system, based on Eqs. ([I]) and ([2|, were performed by varying the nonlinearity 
coefficient, a, and (as above) the coupling coefficient, 7. The initial conditions were taken in the form of Eq. (32), 
unless stated otherwise. 
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FIG. 1: (a) Absolute values \u(z,t)\ and \v(z,t)\, as functions of the propagation distance, z, and temporal 
coordinate, t, obtained from the numerical solution of the linear version of Eqs. ([TJ) and ([ 2 ]) with initial conditions 
( [32| , in the case of weak coupling, 7 = 0.1, q = 1. (b) The evolution of the integral energy of the v component, 
defined as per Eq. 0 and normalized to the initial energy. The evolution of energy of the u component is 

essentially the same, according to Eq. (|5j). 


We start by addressing the weakly coupled system with weak nonlinearity, viz., the one with 0 < d < 1 and 
0 < 7 <C q = 1. For 7 = 0.1 and a = 0.1, Figs. [5^a) and (b) demonstrate that the focusing nonlinearity readily causes 
self-trapping of a robust oscillating quasi-soliton. Thus, the weak nonlinearity, while breaking the CV symmetry (see 
above), creates the self-confined modes. 

Next, we increase the strength of the coupling to 7 = 0.8, keeping the nonlinear term small, with a = 0.1. In 
this case, Figs. [6|a) and (b) demonstrate strong self-focusing of the modes occurring around z = 10, followed by the 
propagation of the confined mode in a sufficiently robust form, although with more conspicuous emission of radiation 
waves than in the case of 7 = 0.1, cf. Fig. [5ja). Thus, in this case too, the system tends to form oscillatory 
quasi-soliton modes. 

The formation of these solitons is readily explained by Eq. (21) for u. Indeed, it is easy to check that the width and 
amplitude of the emerging solitons satisfy conditions (18). The solitons are observed in Figs. [5] and [6] in an oscillatory 
form, which is different from the stationary solution (24), in accordance with the well-knows fact that perturbed NLS 
solitons may feature long-lived vibrations, similar to those observed in these figures mm- 

For 7 > 0.9 and the same weak nonlinearity, <7 = 0.1, the simulations demonstrate that amplitudes of both modes, 
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FIG. 2: The same as in Fig. [lj but for the strong coupling, 7 = 0.9. 


u and v, diverge after a short propagation distance, which implies that the symmetry breaking takes place in these 
cases, that are close to the threshold (11) of the symmetry breaking. The small difference of 7 = 0.9 from the exact 
linear threshold, 7 = 1 , is compensated in this case by the presence of the nonlinearity which, as said above, is also a 
CP-symmetry-breaking factor. 

The swap of initial conditions | 
simulations, performed for input 


and (|33j) in the nonlinear system produces a strong effect. Indeed, in the above 
, the pulse was launched in component it, where the nonlinearity is self-focusing 
[see Eq. 0], while initial conditions ( |33| ) imply that the pulse is launched into component v with the self-defocusing 
cubic term, see Eq. ([ 2 ]). Accordingly, in the latter case, the simulations produce the results displayed in Fig. [ 7 } instead 
of the quick self-trapping (cf. Fig. [6]), the pulse launched in the v component features slow expansion. An additional 
difference is that the frequency of oscillations observed in the latter case is approximately half of that observed in 
Fig. [§ 

The increase of cr at a fixed coupling constant, 7 , enhances the CP-symmetry-breaking effects, and eventually leads 
to destruction of the quasi-soliton. In particular, for the weakly coupled system, with 7 = 0.1, the destabilization of 
the quasi-soliton sets in at critical value u = 0.5, as shown in Fig. [5] In this case, the integral energy slowly grows 
with z, which is followed by blowup at very large values of z (not shown here in detail). 

By means of systematical simulations, we have collected the critical values of cr, at which the quasi-soliton suffers 
the onset of the destabilization, eventually leading to the blowup, at increasing values of the coupling constant, 7 . 
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FIG. 3: The same as in Fig. [5J but for initial conditions (33), with the difference that panel (b) displays the 

evolution of the energy in the u component. 


The corresponding dependence between a and 7 , shown in Fig. (§, naturally demonstrates that the critical strength 
of the nonlinearity vanishes when 7 approaches the threshold of the symmetry breaking in the linear system, 7 = 1 , 
see Eq. (11) [recall the normalization is fixed by setting q = 1 in Eqs. |T]) and Q]. 

The analytical approximation based on Eqs. (19) and (24) was tested for the broad solitons too. Figure 10 


shows that the respective analytical and the numerical results are almost identical, thus validating the analytical 
approximation for the nonlinear system. 


C. Stationary gap solitons in the nonlinear system 


The quasi-solitons considered above are built as breathers, featuring permanent oscillations in both components. 
On the other hand, Eqs. ( 28p0 ) predict the existence of stationary gap solitons in the same system. To check this 
possibility in the numerical form, we solved Eqs. ( 26|27 ) by means of the Newton’s method [52]. This was done using 
the approximate analytical solution, given by Eqs. (28) and (30), as the initial guess. The results are produced here 
for cr = l and different values of the coupling constant, 7. 

In the case of weak coupling, 7 = 0.1, when it is natural to expect the solutions to be strongly asymmetric, in 
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FIG. 4: (a) The numerical simulations for absolute values \u(z,t)\ and k>(£,i)| in the weak-coupling regime, with 
7 = 0.1, q = 1 and the initial conditions taken as per Eqs. |l9| ) and (23) (with the upper sign) at z = 0. (b) The 

respective analytical solutions. 


terms of the two components, the numerical solution at k = 0 , i.e., at the center of the bandgap, is very close to its 
above-mentioned analytical counterpart, as seen in Fig. EH For larger values of 7 , the numerical solution differs from 
the analytical one, obtained under condition 7 <C q = 1, although the difference remains relatively small for 7 = 0.5, 
as shown in Fig. 12 The difference becomes essential for 7 = 0.9 [in fact, very close to the symmetry-breaking 
threshold & as can be seen in Fig. [l3] 

At k ^ 0, the numerically found solutions are still close to the the strongly asymmetric analytical expressions given 
and (30) for k = 0, provided that \k\ is small enough, see Fig. 


by Eqs. ([28]) and (30) for k = 0, provided that |fc| is small enough, see Fig. 14 for 7 = 0.2 and k = 0.2. However, 
at larger k, such as k = 0.8 with the same 7 = 0.2 [note that k = 0.8 falls into the bandgap ( 12 ) in this case], the 
numerical solution for the V component strongly deviates from the analytical expression given by Eq. (27) for k = 0, 
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FIG. 5: The same as in Fig. [l] but produced by simulations of the full nonlinear system § § with q = 1,7 = 0 . 1 , 
and a = 0.1. The weak radiation field around the emerging quasi-soliton is virtually invisible, if local powers, 
\u(z,t)\ 2 and \v(z,t)\ 2 , are displayed instead of the amplitudes, \u(z,t)\ and \v(z,t)\. 


while the U component is still close to the simple analytical form (28 1 , see Fig. 15 


An analytical solution for strongly asymmetric gap solitons (corresponding to 7 <C 1) was also obtained above in 
the form of Eqs. (28) and (31) for k = —(3/5 )q. For 7 = 0.1, this solution is virtually identical to its numerically 
found counterpart, as shown in Fig. |16[ 

Finally, the stability of the numerically generated gap solitons was tested by using them as initial conditions in 
direct simulations of Eqs. 0 and 0. The results, not shown here in detail, demonstrate that all the tested examples 
of gap solitons are stable, both at the center of the gap, k = 0 , and off the center, including values of the coupling 
constant (such as 7 = 0.9) which are close to the symmetry-breaking threshold ©• 


IV. CONCLUSIONS 

The objective of this work is to introduce a model of a dual-core waveguide which may implement an optical system 
featuring the CV symmetry. Essential ingredients of the model are opposite GVD signs in the two cores, a phase- 
velocity mismatch between them, and the linear coupling of the gain-loss type, which is possible if the waveguide is 
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FIG. 6 : The same as in Fig. [5j but in the case of strong coupling, 7 = 0.8. 


embedded into an active medium, or may be provided by the propagation in the medium of the Type-II (three- 
wave) type, neglecting the depletion of the SH (second-harmonic) pump. Nonlinear cubic terms, which destroy the 
symmetry, were considered as well (in the case of the y ^ medium, they are different from those considered here). 
It is predicted in an approximate analytical form and demonstrated numerically that the linear system gives rise 
to expanding Gaussian pulses. Relatively weak nonlinearity produces an essential effect, building broad oscillatory 
quasi-solitons, which are destroyed in direct simulations if the nonlinearity is too strong. Further, the analysis predicts 
a general family of stationary gap solitons in the nonlinear system, that have been also found and checked for the 
stability in the numerical form, the broad solitary pulses being a limit case of the gap solitons near the bottom edge 
of the bandgap. 

The analysis can be continued by considering higher-order modes [it is well known that linear Schrodinger equations 
(21) give rise to higher modes in the form of Hermite-Gauss wave functions] and interactions between fundamental 
solitons in the nonlinear version of the system. On the other hand, it was mentioned above that Eq. (21) for v 
suggests the existence of dark solitons in the present system, which is an interesting issue too. Still another possibility 
is to consider the different form of the cubic nonlinearity, corresponding to the underlying y( 2 ) system. A challenging 
extension is to consider a two-dimensional version of the model, which may be based on a dual-core planar waveguide 
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FIG. 7: The same as in Fig. [6j but for initial conditions (33). 


embedded into the active medium. 
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FIG. 11: Numerical (“Newton” - N) and analytical solutions (A) for the two components of a strongly asymmetric 
gap soliton, U(t) and V(t) (note different scales in the panels), obtained from Eqs. (26) and (27) for q = 1, a = 1, 
7 = 0.1, and k = 0. The respective analytical solution is given by Eq. (28) and (30). 




FIG. 12: The same as in Fig. 11 but for 7 = 0.5. 
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FIG. 14: The same as in Fig. 
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but for 7 = 0 . 2 , and the numerical solution taken for k = 0.2 [recall that analytical 
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[47] G. Agrawal, Nonlinear Fiber Optics, Optics and Photonics (Elsevier Science, 2001). 

[48] E. N. Tsoy, I. M. Allayarov, and F. K. Abdullaev, Opt. Lett. 39, 4215 (2014). 

[49] V. V. Konotop and D. A. Zezyulin, Opt. Lett. 39, 5535 (2014). 

[50] S. Nixon and J. Yang, ArXiv e-prints (2014), 1412.6113. 

[51] C. E. Ruter, K. G. Makris, R. El-Ganainy, D. N. Christodoulides, M. Segev, and D. Kip, 6, 192 (2010). 

[52] G. I. Stegeman, D. J. Hagan, and L. Torner, Opt. Quant. Elect. 28, 1691 (1996), 

[53] and U. Peschel, Progr. Opt. 41, 483 (2000). 

[54] A. V. Buryak, P. Di Trapani, D. V. Skryabin, and S. Trillo, Phys. Rep. 370, 63 (2002). 

[55] H. Suchowski, G. Porat, and A. Arie, Laser Opt. Rev. 8, 333 (1014). 

[56] Y. V. Kartashov, B. A. Malomed, and L. Torner, ArXiv e-prints (2014), 1408.6174. 

[57] C. M. de Sterke and J. E. Sipe, Progress in Optics XXXIII 33, 203 (1994). 

[58] B. J. Eggleton, C. M. de Sterke, and R. E. Slusher, J. Opt. Soc. Am. B 16, 587 (1999). 

[59] R. Blit and B. A. Malomed, Phys. Rev. A 86, 043841 (2012). 

[60] D. Anderson, Phys. Rev. A 27, 3135 (1983). 

[61] B. A. Malomed (Elsevier, 2002), vol. 43 of Progress in Optics, pp. 71 193. 

[62] M. Davis, Numerical methods and modeling for chemical engineers, Wiley series in chemical engineering (Wiley, 1984). 



























18 





FIG. 16: Numerical and analytical solutions for the two components of a strongly asymmetric gap soliton, U(t) and 
V(t), obtained from Eqs. (26) and (27) for q = 1 , er = 1 , 7 = 0.1, and k = —(3/5 )q. The respective analytical 

solution is given by Eqs. (28) and (31). 



























